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Preface 


The  purpose  of  this  study  was  to  experimentally  and 
analytically  determine  certain  aerodynamic  characteristics 
of  a recently  proposed  high  subsonic,  forward  swept  v/ing, 
with  a new  supercritical  airfoil  section.  It  was  my 
intention  to  evaluate  the  effectiveness  of  this  section 
and  to  compare  the  aerodynamic  characteristics  of  a for- 
ward swept  wing  to  a similar  aft  swept  wing.  The  root 
chord,  tip  chord,  wing  span,  and  wing  area  was  held  con- 
stant for  each  wing. 
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Abstract 

This  study  consisted  of  modeling  eind  wind  tunnel 
testing  of  a high  speed,  subsonic,  low  aspect  ratio,  for- 
ward swept  wing  with  an  advanced  supercritical  airfoil 
section  for  the  purpose  of  determining  its  lift,  drag,  and 
pitching  moment  characteristics  as  compared  to  a similar 
aft  swept  wing.  Tests  were  conducted  at  Mach  numbers  of 
0.63  to  0.93  in  the  Air  Force  Flight  Dynamics  Laboratory's 
Trisonic  Gasdynamic  Facility  located  at  Wright-Patterson 

f 

Air  Force  Base,  Ohio.  Two  wing  configurations,  forward 
and  aft  swept,  were  tested  and  compared  to  computer  pre- 
dictions provided  by  the  Unified  Subsonic-Supersonic 
Program  (Woodward's  Version  B).  The  results  indicated 
that  the  forward  swept  wing  was  capable  of  higher  useable 
angles  of  attack  while  maintaining  a lower  drag  coefficient 
for  angles  of  attack  below  eight  degrees.  Wind  tunnel 
test  results  are  presented  in  graphical  and  tabular  form 
for  use  in  future  design  studies  of  similar  aerodynamic 
configurations . 


y AN  AERODYNAMIC. INVESTIGATION  OF  A 
FORWARD  SWEPT  WING 

I , Introduction 


Background 

This  study  addresses  the  problem  of  determining  aero- 
dynamic characteristics  of  a forward  swept  wing  vehicle 
proposed  by  the  Aii  Force  Plight  Dynamics  Laboratory(AFFDL) • 
In  the  past,  the  use  of  a forward  swept  wing  was  generally 
ruled  out  due  to  the  increased  root  loads  requiring  sub- 
stantial increases  in  structural  weight.  With  the  advent 
of  advanced  composites,  it  is  possible  to  negate  these 
weight  penalties  (Ref  A study,  of  a forward  swept 

wing  was  needed  to  obtain  a data  base  for  use  in  the  design 
of  forward  swept  wing  aircraft.  The  proposed  wing  design 
was  generated  to  be  compatible  with  an  existing  half-span 
fuselage  from  AFFDL.  Plan  and  side  view  model  schematics 
are  shown  in  Figs  1 and  2. 

Problem 

The  AFFDL  was  particularly  interested  in  the  aero- 
dynamic investigation  of  a highly  cambered,  supercritical 
airfoil  on  a forward  swept  wing  design  as  compared  to  a 
similar  aft  swept  wing  design.  To  maintain  similar  geom- 
etry on  each  wing,  the  tip  chord,  root  chord,  wing  span, 

and  Wing  area  were  held  constant.  The  following  informa- 
t;i  on  was  required  i 


1.  Determination  of  the  effect  of  cweep  on  the  aero- 
dynamic characteristics. 

2.  Determination  of  the  interference  between  the 
body  and  the  swept  wing. 

The  cambered  wing  was  considered  beneficial  in  improving 
the  stalling  characteristics  of  forward  swept  wings  by 
delaying  leading  edge  separation  (Ref  6 ill).  The  extent 
of  the  anticipated  effect  on  lift  and  drag  was  to  be 
determined. 


Scope 

The  forward  swept  wing  was  tested  and  compared  to  a 
similar  aft  swept  wing  over  numerous  high  subsonic  speeds. 
Each  model  was  analyzed  over  an  extensive  range  of  angles 
of  attack.  The  proposed  model  and  test  parameter  ranges 
included! 


Sweep  at  Quarter  Chord 
Taper  Ratio 
Mach  Number 
Reynolds  Number 
Angle  of  Attack 


-45.0  and  +45.0  (degrees) 
0.4 

0.63  to  0.93 

2,000,000  (per  foot) 
-5.0  to  +15.0  (degrees) 


Approach 

Wind  tunnel  teats  were  the  primary  means  of  determin- 
ing the  lift,  drag,  and  pitching  moment  coefficients  on  each 
model.  The  APFDL’s  two  foot,  Trisonic  Gasdynamic  Facility 
(TGF) , located  at  Wright-Pattcrson  Air  Force  Base,  Ohio,  was 
solocted  for  the  tests.  Data  reduction  was  accomplished 


simultaneously  during  tho  tests#  via  the  Control  Data 
Corporation  (CDC)  6600  digital  computer  tied  directly 

to  the  TGF.  The  force  and  moment  data  reduction  program 

was  provided  by  AFFDL. 

The  Unified  Subsonic -Supersonic  Aerodynamics  Program 
(USS  Aero)  provided  analytical  predictions  which  were 
utilized  in  predicting  the  wind  tvinnol  results.  The 
computer  results  were  used  to  estimate  aerodynamic  char- 
acteristics for  each  model  for  stress  analysis  purposes. 

The  wind  timnel  test  data  would  be  compared  to  the  USS  Aero 
Program  predictions  to  indicate  the  accuracy  of  this  program. 

Materials  and  Equipment 

For  this  study  the  AFFDL  half-span  fuselage  section 
was  used  in  conjunction  with  two  swept  wings.  Each  wing 
was  mounted  at  the  same  root  chord  position  on  the  fuse- 
lage and  constructed  using  a 5*5  percent  thick,  cambered, 
supercritical  airfoil  provided  by  AFFDL.  The  geometry  of 
each  wing  was  kept  as  similar  as  possible  to  allow  better 
comparison.  Other  modifications  to  the  existing  modal 
included I 

1.  New  center  body  for  wing  support. 

2.  New  center  fuselage  section  to  allow  fnr 
pressure  linos. 

3.  Thirteen  pressure  taps  in  each  wing. 
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II . Facility  And  Model  Deocription 


Facility  Used  for  Teatinfi 

The  Air  Force  Plight  Dynamics  laboratory 's  Trisonio 
Gaadynamic  Facility  (TOP)  was  utilized  for  the  wind  tunnel 
tests.  This  facility  was  chosen  for  several  reasons • The 
TOP  is  a closed  circuit,  variable  density  wind  tunnel 
capable  of  continuous  flow.  It  may  bo  operated  at  sub- 
sonic, transonic,  and  supersonic  speeds  through  a range  of 
<!• 

Mach  numbers  from  O.23  to  4.76.  The  facility's  axia],  flow 
compressor  was  maintained  at  constant  RPM  with  the  use  of 
two  synchronized  motors,  capable  of  up  to  11,500  horse- 
power, Hydraulically  actuated  stator  blades  controlled, 
instantaneously , the  volume  flow  and  pressure  rise,  Water  ■ 
cooled  heat  exchangers  maintained  a stagnation  temperature 
of  100  P + 1 F, while  two  dryer  systems  maintained  a dew- 
point necessary  to  prevent  llquofacation  of  water  vapor  in 
the  test  section,  Quantitative  data  on  test  section  turbu- 
lence levelp  wrrti  j.  ri,  available,  but  a honeycomb  and  screen 
arrangement  located  in  the  tunnel  stagnation  section  were 
used  to  reduce  turbulence  in  the  test  section.  The  variable 
density  capability  allowed  for  a constant  high  unit  Reynolds 
Number,  approximately  2.0  million  per  foot,  to  be  used  as 
the  Mach  number  varied.  The  charactoristics  of  the  tran- 
sonic test  section  are  given  in  Table  I . 
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T.ablQ  I 

TGF  Transonic  Teat  Section  Characteristics 


I 


- 1 

'f- 

■I 


' i 


I 


i 


Mach  Number  Range 
<l  PSF  (Max) 

Pq  PSF  (Max) 

Pq  PSF  (Max) 

R/Ft  (Max) 

To  (F) 

Cross  Section  (in.) 
Density  Altitude  (ft.) 
Number  of  Slots 


0.30-1,20 

1,480 

4,000 

2,110 

8,250,000 

100 

15  X 15 
50,000 
56 

(Ref  12i3) 


An  11.426  inch  diameter  window  in  the  test  section 
permitted  access  to  the  model  for  rapid  configuration 
changes  an  well  ae  an  unobstructed  view  of  the  model  during 
testing,  The  window  facilitated  the  use  of  oil  flow  photo- 
graphy as  a flow  visualizatdv^n  method, 

The  model  was  tupported  by  the  Uenora.i  Dyn-'rA 
span  rig  mounted  with  the  center  of  rotation  in  the  center 

4 

of  the  test  section  viewing  window.  The  pitch  range  was 
negative  five  to  positive  fifteen  degrees  wi'^h  on  accuracy 
of  + 0.01  degrees.  While  maintaining  the  desired  Mach 
number,  the  wide  range  of  possible  angles  of  attack  was 
varied  (Ref  12il-17) . 

The  model  support  system  allowed  the  use  of  a five 


4 

.} 


i 

» 

j 


i 

i 


I 


I 


I J' 


I 


i s* 
( i 

1; 


I 


7 


component  internal  typo  strain  gage  balance  to  obtain 
force  and  moment  data.  The  balance  oelcctod  for  testing 
was  a General  Dynamic  half-span  balance,  number  2M-259. 

It  was  capable  of  measuring  normal  and  axial  forces  of 
250  and  30  pounds  respectively.  It  allowed  for  roll,  yaw, 
and  pitch  moments  of  500  inch-pounds  to  be  determined 
(Ref  3) . The  load  sensing  elements  of  the  balance  were 
connected  to  a Wheatstone  bridge  arranged  in  such  a manner 
as  to  give  an  electrical  signal  proportional  to  the  defor- 
mation. The  electrical  signal  in  millivolts  was  fed  to  self- 
balancing  potentiometers  and  by  use  of  shaft-to-anolog-to- 
dlgital ' converters  were  translated  to  digital  outputs . The 
outputs  were  reduced  to  a more  usable  aerodynamic  coeffi- 
cient form  by  a CDC  6600  computer.  These  coefficients 
were  tabulated  and  plotted.  This  force  measurement  system, 
therefore,  provided  Instantaneous  visual  monitoring  and  an 
automatic  digital  recording  capability  (Ref  12«10-17). 

The  pressure  sensing  instrumentation  consisted  of  differ- 
ential pressure  transducers.  Each  transducer  was  capable 
of  reading  up  to  48  separate  pressures  as  each  tap  was 
electrically  stlooted,  A series  of  calibration  tests 
Indicated  the  balance  and  associated  instrumentation  cap- 
able of  measuring  forces  and  moments  within  + 1 porcont  of 
the  applied  Ick  d.  The  pressure  data  was  calibrated  to 
wli.hin  one  tenth  of  a pound  per  square  foot. 


a 


Modolfl  Tented 

For  this  series  of  tests,  two  entirely  now  wings  were 
constructed  for  use  with  existing  General  Dynamics  half-span 
fuselage.  The  geometry  requirements  for  each  wing  were 
dictated  by  the  use  of  the  half-span  fuselage  due  to  con- 
straints in  available  construction  techniques  and  time 
available  I A new  fuselage  center  section  was  constructed 
to  allow  support  of  each  wing  at  the  same  fuselage  station 
thus  requiring  minimum  time  for  the  interchange  of  the 
wings  during  testing.  As  a result,  the  forward  swept  con- 
figuration tested  had  less  than  ideal  pitching  moment 
characteristics.  Pertinent  model  dimensions  are  given  in 
Table  2 and  the  models  are  shown  in  Figs  3 and  4, 

Table  II 
Model  Geometry 


Wing  Planform  Area  (sq.  in.) 

Root  Chord  (in.) 

Tip  Chord  (in.) 

Taper  Ratio 

Wing  Semi -Span  (in.) 

Sweep  Angle  of  Quarter 
Chord  (degrees) 

Aspect  Ratio 

Mean  Aerodynamic  Chord  (in.) 


25.725 

6.125 

2.45 

0.40 
7. '58 

+45,  -45 
3.34 
5.62 


The  airfoil  to  bo  utilized  was  a nowly  developed, high 


camber  I supercritical  airfoil  section  of  5*5  percent  local 
chord  thickness  ratio.  I'ho  maximum  chord  of  6.125  inches 
occurred  at  the  root  thus  making  the  thickness  of  the 
section  0,337  Inches.  For  support  purposes,  an  eighth  inch 
thick  steel  plate  was  embedded  in  an  epoxy  mold  of  each 
wing  and  subsequently  mounted  to  the  fuselage  center  sec- 
tion. The  chord  on  'each  wing  was  aligned  with  the  flow 
direction  thus  allowing  support  of  the  wings  on  the  same 
fuselage  center  section.  This'  airfoil  was  used  due  to  its 
high  lift  properties  in  transonic  flow  and  for  the  reduction 
of  leading  edge  separation  it  exhibited.  By  utilizing 
identical  sectional  and  planform  geometry,  the  data  com- 
parison and  model  design  and  construction  of  the  two  wings 
could  be  greatly  accelerated.  Although  several  significant 
changes  were  made  to  the  final  design  of  the  models,  the 
basic  half-span  fuselage  with  available  ZM-259  balance 
could  be  used  with  the  newly  constructed  wings  (Ref  5) • 
Pressures  were  taken  at  thirteen  separate  locations 
on  eeoh  wing  in  one  chordwise  cut.  Tho  pressure  taps  were 
aligned  at  a mid-span  location  to  limit  the  amount  of  both 
fuselage  and  tip  Interference  on  the  pressure  data,  in 
Table  3»  these  pressure  tap  locations  are  listed  in 
percent  chord  for  the  4.3  inch  local  chord  position  on 
each  wing.  Only  taps  2 and  3 are  located  on  the  lower 
surface  of  the  airfoil. 


Table  III 


Pressure  Tap  Locations 


r ■ 

1 Tap  Number 

r 

fc  Chord 

1 • 
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1 ' 

.v 

1 i 1 
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Experimental  Procedures 


Force  Measurements 

The  reduction  of  data  using  the  half-span  balance 
output  voltage  was  performed  on  a CDC  6600  digital  com- 
puter (Ref  12il0).  Aerodynamic  characteristics  are  pre- 
sented in  the  form  of  coefficients  in  the  wind  axis 
system.  The  longitudinal  moments  were  computed  about 
the  balance  center  to  enable  direct  comparison  of  the 
different  wing  configurations.  The  mean  aerodynamic 
chord  of  5*62  inches  was  computed  usingi 

C * I . 

where  Cj.  is  the  root  chord  and  >\  is  the  taper  ratio. 

It  was  used  as  a reference  length  in  the  data  reduction 
program. 

An  important  consideration  in  obtaining  true  axial 
force  measurements  in  most  wind  tunnel  tests  is  the 
treatment  :f  base  drag.  This  drag  arises  from  the  pres- 
sure acting  on  the  blvint  base  of  the  model.  The  base 
pressure  is  influenced  considerably  by  the  presence  of 
the  tunnel  wall,  acting  as  a boundary  to  the  model  base. 
This  influence  normally  causes  the  drag  values  to  be 
smaller  than  they  would  be  if  the  model  could  \>e  tested 
without  a boundary  and  balance  in  place  (Ref  8i4l3). 
Therefore,  the  measurements  taken  in  testing  each  model 
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were  corrected  to  remove  this  base  pressure  effect  as 
followsi 


F.C  - F.  - StC  P-PJ 

where  Fac  is  the  axial  force  adjusted  to  zero  base  drag# 
Fa  is  the  axial  force  measured#  is  the  model  base 
area#  P is  the  free  stream  static  pressure,  and  P^,  is  the 
model  base  pressure  (Ref  9«323)*  The  base  pressure  was 
measured  with  a pressure  probe  located  in  the  base  of 
the  model . 

Test  Procedures 

The  wind  tunnel  tests  were  divided  into  two  main 
procedures.  The  first  test  was  to  obtain  the  aerodynamic 
characteristics  of  each  configuration  to  include  lift# 
drag,  and  pitching  moment  coefficients.  The  method  of 
obtaining  these  characteristics  was  to  vary  angle  of 
attack  over  a range  of  negative  four  to  positive  fifteen 
degrees  while  maintaining  a constant  Mach  number.  Each 
model  was  placed  in  position  with  tunnel  off  to  obtain 
static  pressure  and  zero  tare  readings.  Once  the  tunnel 
was  stabilized  at  the  desired  Mach  number',  these  readings 
were  again  taken  and  corrected  for  wind  off  conditions. 
Numerous  Mach  numbers  over  a range  of  0.63  to  0.93  were 
tested.  The  total  pressure  was  varied  with  Mach  number 
to  maintain  a unit  Reynolds  number  of  about  two  million. 
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The  test  settings  are  given  in  Table  4. 

After  force  measurements  on  both  configurations 
were  completedi  a brief  flow  visualization  study  was  made. 

A white  tinted  oil  was  brushed  on  the  models  for  this 
purpose.  Flight  conditions  representing  those  tested  were 
then  investigat  1.  The  oil  flow  permitted  visualization  of 
the  interactions  of  s'treamlines.  Regions  of  flow  separa- 
tion are  normally  detechable  using  this  method  of  flow 
visualization. 


Table  IV. 
Test  Conditions 


Tq  s 100  P 

Data  taken  at  each  Mach  number  for  alpha  range  of 
-5.0  to  +15.0  degrees 

Pq  varied  to  maintain  R = 2,000,000  per  ft. 
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IV.  Analytical  Procedures 

In  order  to  obtain  initial  predictions  of  the  wind 
tunnel  datai  a new  method  of  calculating  the  pressure 
distribution  and  aerodynamic  characteristics  of  wing- 
body  -tail  combinations  was  used.  The  Unified  Subsonic- 
Supersonic  Aerodynamics  Program  (USS  Aero)  performed 
the  numerical  calculations.  This  program  uses  subsonic 
potential  flow  theory  (Ref  13i2). 

Both  the  forward  and  aft  swept  configurations  were 

modeled  by  subdividing  the  wing  as  well  as  the  fuselage 

into  a large  number  of  panels.  Each  of.  these  panels 

contained  an  aerodynamic  singularity  distribution.  The 

modeling  for  the  aft  and  forward  swept  configurations  is 

shown  in  Figs  5 and  6 » respectively.  A constant  source 

distribution  was  used  on  the  body  panels  and  a vortex 

distribution  having  a linear  variation  in  the  streamwise 

direction  was  used  on  the  wing  panels.  The  pleuiar  option, 

which  uses  control  points  in  the  plane  of  the  wing,  was 

used  and  yielded  accurate  comparisons  with  experimental- 

* 

data.  The  non-planar  option,  which  uses  control  points 
on  the  surface  of  the  wing,  required  double  the  number  of 
panels  and  yielded  only  slightly  better  results.  This 
mode  was  not  used  due  to  'he  large  increase  in  computer 
time  with  minimal  increase  in  accuracy.  The  normal  comp- 
onents of  the  velocity  Induced  at  specified  control  points 


I 'i 
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Fifjure  6.  Forward  Swopt  Model  Used  For  LiCS  Aero  Prediction 
Pro{fram . 
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by  each  cingularlty  distribution  were  calculated  and  make 
up  the  coefficients  of  a system  of  linear  equations  re- 
lating the  strengths  of  the  singularities  to  the  magnitude 
of  the  normal  velocities. 

The  singularity  strengths  which  satisfy  the  boundary 
condition  of  tangential  flow  ‘at  the  control  points  for  a 

■given -Mach- jaiuabftff...and  of  attack  were  determined  by 

solving  the  system  of  equations  using  8.]ir  I'ferat’l'vtJ-pre- 

cedure.  Once  the  singularity  strengths  are  known,  the 
pressure  coefficients  were  calculated,  and  the  forces  and 
moments  acting  on  each  oonf^ "uxation  were  determined  by 
numerical  integration  (Rt  ) . 

A skin  friction  coeffig.^ant , Cf,  had  to  be  determined 
and  subtracted  from  the  experimental  drag  coefficient  to 
obtain  a valid  comparison  of  the  relative  accuracy  of 
Woodward's  USS  Aero  program.  This  allowed  comparison  of 
the  induced  drag  coefficient  computed  by  the  potential 
flow  analysis  to  the  experimental  drag  coefficient  without 
skin  friction. 

The  skin  friction  drag  of  each  configuration  was  cal- 
culated from  the  T'  method  (Ref  7il2).  The  T'  method  ds 
based  on  the  calculation  of  a compresslbie  skin  friction 
coefficient  from  a reference  skin  friction  coefficient,  for 
a selected  Mach  number,  Reynolds  number,  and  adiabatic 
wall  temperature,  Smooth  flat  plate , adiabatic . wall , and 
turbulent  boundary  layer  conditions  were  assumed.  In 
addition,  transition  was  assumed  to  occur  at  tho  loading 
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edge  of  each  configuration  component.  Most  of  the  program 
is  involved  with  computing  wetted  ar^^as  and  rofex’ence 
lengths  for  each  component.  Since  both  configurations  had 
identical  wetted  areas,  Table  5 represents  the  skin  fric- 
tion coefficient  calculated  for  each  Mach  number  tested. 
These  values  of  skin  friction  were  subtracted  from  experi- 
mental data  before  comparison  was  made  to  theoretical  data. 


Table  V. 

Skin  Friction  Coefficients  For 
Test  Mach  Numbers 


Mach 

Skin  Friction 

Coefficient 

0.63 

0.01532 

0,70 

0.01523 

0.80 

0.01509 

0.85 

0.01502 

0.90 

O.OlW 

0.93 

0.01489 

Three  Mach  numbers,  representative  'of  the  test  Mach 

I 

numbers,  were  utilized  to  obtain  theoretical  data  for 
comparison  to  the  tunnel  data  for  each  configuration.  Good 
correlation  between  theory  and  experimental  data  was 
achieved . 
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V.  Reaulta 

The  wind  tunnel  investigation  of  the  forward  and 
aft  swept  wing  models  resulted  in  a series  of  aerodynamic 
ooeffloientB  for  each  configuration  tested  • The  coefficients 
determined  for  each  model  included  lift,  drag,  and  pitching 
moment  coefficients , The  coefficients  are  presented  in 
the  wind  axis  system  as  shown  in  Fig  ?•  The  results 
discussed  in  this  section  are  for  Mach  numbers  of  0i63 
and  0.93.  The  results  obtained  for  Mach  numbers  between 
these  two  values  are  presented  in  graphical  form  in 
Appendix  B,  In  addition,  flow  visualization  photographs  at 
flight  Conditions  typical  of  those  tested  are  included, 

Uft 

The  coefficients  for  the  forward  and  aft  swept  models 
are  plotted  as  a function  of  angle  of  attack  in  Pigs  8 and 
9 for  the  low  and  high  Mach  number  tests,  respectively. 

These  values  are  typical  of  those  obtained  at  all  remaining 


Mach  numbers  tested.  The  slope  of  the  lift  curve  is 
approimately  0.060  to  0.O65  per  degree,  depending  on  the 
angle  of  attack  and  Mach  number  tested.  For  all  Mach 
numbers  tested,  the  aft  swept  wing  produced  a higher  lift 
coefficient  for  a given  angle  of  attack  as  compared  to  the 


forward  swept  wing.  As  the  Mach  number  was  Increased  to 
high  subsonic  values,  the  aft  swept  model  lift  curves 
tended  to  shift  toward  the  forward  swept  model  lift  curves 
until  the  two  curves  varied  by  ^ equal  to  0.035  at 


PlGure  7.  Wind  Axiei  Coordinate  uood  for 
Aorodynamic  Coefficients 


MflCH  =0.93 


I 


i 


l!  ! 

' i 

^1  f 


M s=  0.93.  The  chance  in  angle  of  aero  lift  was  the  major 
contributor  to  the  shift  of  the  aft  swept  model  lift  curve. 

At  M = 0.63  the  zero  lift  angle  for  the  aft  swept  model  was 
-3.2  degrees  whereas  at  M = O.93,  it  became  -2.0  degrees. 

The  forward  swept  model  had  an  angle  of  zero  lift  which 
varied  between  -2.0  and  -1.6  degrees.  The  aft  swept  .aodel 
showed  signs  of  stalling  at  a lower  angle  of  attack  than 

1 

the  forward  swept  model  for  all  Mach  numbers  tested. 

j 

I 

Drag 

A typical  set  of  drag  coefficients  is  plotted  in 
Figs  10  and  11  as  a function  of  angle  of  attack  and  in  ; 

Figs  12  and  13  as  a function  of  lift  cpefficient.  For  M = 1 

0.63,  the  drag  coefficient  at  zero  angle  of  attack,  Cdq  » 
was  0.038  for  the  forward  swept  model  and  0.042  for  the  | 

aft  swept  model.  As  Mach  number  was  increased  to  M - 0.93»  j 

the  values  of  Cqq  changed  to  0.042  and  0.051  for  the  for- 
ward and  aft  models,  respectively.  For  all  Mach  numbers 

■ i 

tested  the  forward  swept  model  displayed  less  drag  up  to  i 

1 

approximately  ten  degrees  alpha.  Above  this  angle  the  two 
drag  curves  crossed  and  the  forward  swept  model  produced 
more  drag  for  a given  alpha. 

The  maximum  lift  to  drag  ratio  may  be  determined  by 
drawing  the  tangent  to  the  Cd  versus  Cj;^  curve  from  the 
origin  (Ref  9«330).  For  M = O.63,  this  yields,  an  L/D I max 
for  the  aft  swept  model  of  8.3  £^t  a Cl  of  0.51'  For  the 
forward  swept  model,  this  ratio  was  8.4  at  a of  0.40. 
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Fii^ure  10.  Cq  vs  for  M " 0.63 
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This  correoponda  to  an  angle  of  attack  of  4.1  degrees  and 
6.2  degrees  for  the  aft  and  forward  swept  models i respec- 
tively. This  shows  thatvth©  . forward  swept,  model  was  capable 
of  higher  angles  of  attack  for  thie  same  lift  to  drag  ratio. 
The  lift  to  drag  ratio  is  plotted  in  Pigs  l4  and  15  as  a 
function  of  angle  of  attack  for  direct  comparison. 

Figure  16  shows  the  effect  of  Mach  nianber  on  the 
drag  coefficient  at  a constant  lift  coefficient.  A point 
on  the  curve  that  is  sometimes  used  for  reference  is  the 
drag  divergence  Mach  number,  MpD  • The  is  the  Mach 
number  at  which  dC^/dM  is.  equal  to  0.10  (Ref  2i2).  The 
Mdd  decreases  with  an  increase  in  alpha  or  an  increase  in 
lift  coefficient.  For  the  aft  swept  model,  Mdd  is  above 
0,93  for  alpha  less  than  or  equal  to  five  degrees.  It  is 
approximately  0.84  at  an  alpha  equal  to  eight  degrees.  For 
the  forward  swept  model  a drag  divergence  Mach  was  not 
indicated  since  the  drag  continued  to  increase  in  an  approx 
imately  linear  manner  with  increase  in  Mach  number. 

Pitching  Moment 

The  pitching  moment  about  the  balance  center  is 
represented  in  Figs  17  and  18.  The  slope  of  this  curve 
varies  little  with  Mach  number  for  each  model.  Due  to 
constraints  in  model  fabrication  , the  forward  swept  con- 
figuration resulted  in  a positive  slope  to  the  pitching 
moment  curve.  This  indicates  the  need  of  shifting  the 
location  of  the  root  chord  to  a more  rearward  position  to 
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Figure  l4,  L/D  vs  for  M 
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obtain  a atable  pitching  moment  for  tliia  configuration. 
Although  the  slopoa  of  thr>  two  curves  are  opposite  .a  oom- 
parison  of  the  relative  nagnltuda  of  each  slope  indicated 
t'  a f rwurd  swLpu  model  was  less  sensitive  to  changes  in 
angle  of  attack.  For  M •»  0.63,  at  an  angle  of  attack  of 
four  degrees,  the  slope  of  the  pitching  moment  curve  is 
approximately  +0.020  per  degree  for  the  forward  swept  model 
and  -0.038  per  degree  for  the  aft  swept  model,  The  slope 
changes  are.  of  course . dependent  on  the  angle  of  attack. 


The  lift  and  drag  coefficient  curves  for  both  the 
USS  Aero  predictions  and  experimental  data  are  shown  for 
oomparison  in  Pigs  19  to  21  for  the  fofward  swept  model 
and  in  Figs  25  to  2?  for  the  aft  swept  model.  The  lift 
coefficient  curves  compare  extremely  well  with  the  USS 
Aero  program  predicting  approximately  the  same  slope  and 
zero  lift  angle  of  attack  as  obtained  from  experimental 
data.  Drag  comparisons  are  shown  in  Figs  22  to  24  for 
the  forward  swept  model  and  in  Pigs  28  to  30  for  the  aft 
swept  model.  Note  that  the  USS  Aero  program  only  predicts 
induced  drag  and  does  not  include  any  skin  friction  drag. 

I 

The  experimental  data  was  corrected  by  subtracting  out  the 
computed  skin  friction  drag  as  discussed  previously.  The 
predicted  drag  curves  displayed  the  same  shape  of  the  drag 
polar  obtained  from  experimental  values  with  small  diff- 
orencoB  in  values  of  Cdq.  For  the  panel  distribution  used 
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on  these  models i the  USS  Aoro  program  provided  woll- 
bohavod  and  consiatont  data  for  comparison  to  experimental 
results. 

Pressure  Distribution 

The  pressure  data  obtained  from  this  test  are  pre- 
sented as  a pressure  poeffloient,Cp  , plotted  as  a func- 
tion of  percent  chord  as  shovm  in  Figs  31  to  3^*  At  M « 
0.63 1 the  data  produced  a smooth  curve  with  no  shook  pre- 
sent. For  M « 0.93 1 for  the  aft  swept  model,  a definite 
shook  oan  be  observed  between  70  and  75  percent  chord 
location.  For  the  forward  swept  model  a oompresslon  begins 
at  the  70  percent  chord  location  but  is  gradual  and  there- 
fore should  not  be  referred  to  as  a shook. 

Although  not  plotted,  the  five  pressure  taps  located 
around  the  leading  edge  provided  data  which  indicated  that 
even  at  the  highest  angle  of  attack  tested,  leading  edge 
separation  did  not  ooour.  Only  at  negative  alphas  did 
separation  develop  on  the  lower  surface  due  to  the  largo 
amount  of  oamber  in  the  airfoil  tested.  Due  to  the  thinnoss 
of  the  airfoil  section,  more  pressure  taps  on  each  wing  was 
not  possible.  Therefore,  a complete  pressure  distribution 
including  data  between  the  10  and  50  percent  chord  loca- 
tions ware  not  obtained. 

Oil  Visualization 

The  results  of  the  oil  flow  study  can  be  seen  in  Figs 
35  through  40 . Figure  35  shows  the  uft  swept  wing  at 
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<A=  4,  M = 0.63,  Compare  this  to  Fig  36,  the  forward 
swept  wing  at  4,  M - 0.63i  and  it  is  apparent  that 
neither  wing  displays  any  flow  separation.  Some  spanwise 
flow  can  he  seen  on  the  forward  swept  oil  photograph  hut 
it  is  not  known  whether  the  aft  swept  model  produced  span- 
wise  flow  also.  Note  that  a thin  section  on  eaoh  wing 
appearing  about  mid-span  is  due  to  tape  used  to  protect  the 
pressure  taps  and  is  not  an  aerodynamic  phenomena.  Figure 
37  and  Pig  38  present  the  aft  swept  model  at  «5C=  10.  M =• 
0.63  and  the  forward  swept  model  at  the  same  conditions. 

The  wing  tip  of  the  aft  swept  wing  and  the  root  of  the 
forward-  swept  wing  show  indications  of  flow  separation. 

Note  that  the  flow  inward  along  the  trailing  edge  of  the 
forward  swept  wing  is  more  pronounced  than  before. 

Figures  39  and  40  show  the  same  wings  at  10, 

M = 0.93*  The  flow  separation  has  developed  into  a full 
stall  on  the 'tip  of  the  aft  swept  wing  and  on  the  root  of 
the  forward  swept  wing  with  flow  reversal  apparent.  The 
flow  on  the  outboard  section  of  the  forward  swept  wing 
remains  attached  even  at  an  angle  of  ten  degrees. 

The  interference  of  the  body  on  each  wing  could  not 
be  determined  since  photographs  of  each  model  had  to  be 
taken  at  an  angle  from  the  tip.  Plow  patterns  near  the 
root  chord  were  not  clearly  visible  on  cither  model. 
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Figure  37.  Oil  Flov/  of  Forv/ard  Swept  V/ing  at 
M = 0.63  , 10,  Flow  Right  to  Left  • 


Figure  38.  Oil  Flov/  of  Aft  Sv/cpt  Vllng  at 
M 0.63  , 10.  Flow  Right  l;o  Loft 
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VI . Conclusions  and  Recommendations 


Conclusions 

The  following  conclusions  were  drawn  baccd  upon  the 
results  of  the  wind  tunnel  tests  as  discussed  in  Section  V. 

1 . It  was  concluded  that  the  data  obtained  in  this 
study  was  reliable  and  accurate  as  indicated  by  the 
numerous  tests  for  repeatability  and  data  comparisons  with 
the  USS  Aero  Program • 

2.  It  was  concluded  that  the  highly  cambered,  super- 
critical airfoil  alleviated  any  problem  of  leading  edge 
separation  for  this  configuration.  Since  both  models 

did  not  experience  any  leading  edge  separation  effects, 
the  stall  Indications  were  delayed  to  the  point  of  neither 
wing  stalling  even  at  an  angle  of  attack  of  fifteen  degrees. 

3*  It  was  concluded  that  the  forward  swept  wing  was 
capable  of  a higher  useable  angle  of  attack  than  the  aft 
swept  wing.  For  each  Mach  number  tested,  the  lift  coef- 
ficient obtained  for  the  aft  swept  wing  was  matched  by  the 
forward  swept  wing  but,  required  a higher  angle  of  attack. 
When  the  aft  swept  wing  began  to  stall ,,  the  forward  swept 
wing  lift  curve  continued  to  increase  up  to  the  highest 
angle  of  attack  tested.  The  actual  stall  of  each  wing  was 
above  fifteen  degreed  so  the  extent  of  increase  in  useable 
a-nglo  of  attack  could  not  be  determined. 

It  was  concluded  that  the  drag  coefficient  for  the 
forward  swept  wing  was  consistently  lower  than  the  similar 
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aft  swept  wing  for  lift  coefficients  below  approximately 
0,50  which  corresponds  to  an  angle  of  attack  less  than 
eight  degrees  alpha.  This  indicates  that  the  forward 
swept  wing  could  be  used  more  efficiently  than  the  aft 
swept  vdng  for  normal  .cruise  flight  operational  thus 
requiring  less  fuel. 

Recommendations 

The  following  recommendations  follow  from  the  results 
of  the  wind  tunnel  tests 1 

1 . The  USS  Aero  Program  correctly  predicts  the  trends 

in  the  aerodynamic  characteristics  of  the  configurations 

tested.  It  .'is  recommended  that  the  USS  Aero  Program  be 

* 

used  to  predict  the  effects  of  small  design  changes  and 
these  predictions  be  used  to  adjust  the  wind  tunnel  data. 

2.  It  is  recommended  that  the  location  of  the  for- 
ward swept  wing  be  moved  aft  of  the  present  location  to 
obtain  a stable  pitching  moment  configuration.  The  long- 
itudinal stability  coefficients  could  then  be  studied  to 
determine  the  size,  shape,  and  location  of  a canard 
surface . 

3.  The  oil  flow  study  indicated  that  a large  portion 
of  the  root  section  of  the  forward  swept  wing  had  lost  its 
effectiveness  due  to  flow  separation.  It  is  recommended 
that  the  use  of  a forward  swept  canard  in  conjunction  with 
the  forward  swept  wing  be  tested  to  determine  its  effect 
upon  maintaining  attached  flow  at  higher  angles  of  attack. 
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The  use  of  winglets  should  -also  be  considered  to  determine 
their  effect  on  the  aorodynamic  characteristics. 

4.  It  is  recommended  that  the  data  obtained  from 
this  study  be  used  as  a basis  for  a performance  study  of 
a forward  swept  wing  aircraft.  Advantages  and/or  dis- 
advantages of  the  forward  swept  wing  could  be  identified 
to  support  the  development  of  a forward  swept  wing  aircraft. 


iammiiSSSOSSSSS 


1 . Armed  Forces  Supply  Support  Center.  Strength  of 
Metal  Aircraft  Elements . Military  Handbook  5Ai 
March  1965* 

2.  Corning  G.  Supersonic  and  Subsonic  Airplane  Design. 
Ann  Arbor,  Michigani  Braun --Brumfield , Inc.,  1970. 

3.  Klein,  O.J.  Development  of  a Half -Span  Model  Test 

System  for  the  AFFDL-TGF .Son^raot  5-7^-C^3^5  • 

Vol  1 , Ft.  Worth,  Texas I General  Dynamics  Corporation 

Sept.  1975. 

^1-.  Krone,  N.J.  Divergence  Elimination  With  Advanced 
Composites.  AIAA  Paper  75-1009,  August,  1975. 

5.  Lawrence, ' J.R.  Development  of  a Half-Span  Model 
Test  System  For  The  AFFDL-TGF.  Contract  F-33^15-7^- 
C-3085,  Vol.  2,  Ft.  Worth,  Texas 1 General  Dynamics 
Corporation,  September  1975. 

6,  MoCormaokiG.M.  and  W.L.  Cook,  A Study  of  Stall  Phen- 
omena on  a ^ Degrees  Swept  Forward  Wing,  NACA  Tech- 
nical Note  1797.  Washingtoni  National  Advisory 


nical  Note  1797 i Washingtoni  National  Advisory 
Committee  for  Aeronautics,  January  19^9* 

7.  Middleton,  M.D.  and  J.L.  Lundry,  A Computational 
System  For  Aerodynamic  Design  and  Analysis,  NASA 
Contractor  Report  2715,  WashingtonT  National 
Aeronautics  and  Space  Administration,  July  I976. 

8.  Pope,  A.  Wind  Tunnel  Testing.  New  Yorki  John 
Wiley  and  Sons,  Inc.,  1961. 

9.  , and  K.L.  Goin,  High  Speed  Wind  Tumel  Testing. 

New  Yorki  John  Wiley  and  Sons,  Inc.,  1^6$'. 

1 0 . Roskam , J , Flight  Dynamics  of  Rigid  and  Elastic  - 
Airplanes.  Kansas  1 Roskam  Aviation  and  Engineering 
Corporation,  1975. 

11.  Roark,  R.J.  and  W.C.  Young.  Formulas  for  Stress 
and  Strain  (5th  Edition).  New  Yorki  McGraw-Hill 
Book  Company,  1975. 


12.  White,  H.L.,  Triaonic  Gasdynamlc  Facility  User 
Manual.  Task  Niunber  14760^01 . Wright  -Pa  tt  o r son 
Air  Sorce  Base,  Ohloi  Air  Force  Flight  Dynamics 
Laboratory,  June  1973 • 

13.  Woodward,  F.A.  An  Improved  Method  For  the  Aerodynamic 
Analysis  of  Winfi-Bodv-Tail  Oonfigurajiions  in 

tubsonio  and  Supersonic  M.o^  NASA  Contractor  Report- 
228,  Part  i , Waahingtoni  National  Aeronautics  and 
Space  Administration,  May  1973. 


Appendix  A 

Model  Stress  Analysis 


A stress  analysis  of  the  two  newly  constructed  wings 
was  required  to  Insure  they  could  withstand  the  loads 
expected  during  the  wind  tunnel  tests.  A rather  stringent 
ultimate  factor  of  safety  of  4.0  was  required  to  insure 
against  model  failure.  Each  model  was  analyzed  and  eval- 
uated against  this  criterion.  0)wo  basic  assumptions • 
which  tend  to  yield  conservative  estimates,  were  used 
throughout  this  analysis.  First,  the  pressure  distribution 
across  the  wing  was  assumed  uniform.  S.econdly,  the  entire 
load  of  the  wing  was  assumed  concentrated  at  the  centroid 
of  the  wing  instead  of  the  theoretical  center  of  pressure. 
Both  assumptions  produce  higher  bending  moments  and  tip 
loads  than  would  be  expected. 

Design  lLoad,a 

The  maximum  aerodynamic  loads  expected  for  each  wing 

was  determined  using  a value  of  Ct,  =1.0,  which  was 

■^max 

obtained  from  previous  two-dimensional  tests  of  the  super- 
critical airfoil.  The  maximum  aerodynamic  lift  isi 

L = 9 S (1) 

where  q is  the  dynamic  pressure  and  the  reference  area, 

2 

S = 0.2393  ft  . The  maximum  dynamic  pressure  in  this  test, 
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which  occurod  at  M ==  0.93,  is  given  byi 


K P 
z. 


(2) 


where  K is  the  ratio  of  specific  heat,  M is  the  Mach  number, 
and  P is  the  static  pressure.  The  relationship  for  static 
pressure  is  given  byi 

P = (?  C t + liM'  (3) 

5 ^ 

where  Pq  is  the  total  pressure.  Using  Pig  of  Ref  12, 
the  total  pressure  of  the  tunnel  can  be  found  as  a function 
of  Mach  number  and  Reynolds  number.  The  maximum  load 
for  this  test  occurred  at  M = 0.93  and  R = 2,000,000  per 
foot  at  which  Pg  = lOOOpsf.  Using  Eq  (2)  and  Eq  (3),  the 
maximum  dynamic  pressure  was  computed  as  3^6.^  psf. 

Referring  to  Fig  4-1,  the  maximum  aerodynamic  load  is 
shown  to  be  a function  of  angle  of  attack  given  byi 


N “ L 

CoS  <K. 


(4) 


where  N is  the  maximum  normal  aerodynamic  force,  L is  82.9 
Ibf  from  Eq  (1),  and  , the  angle  of  attack, is  fifteen 
degrees.  This  yields  a maximum  aerodynamic  load  of  85.8 


Ibf. 


Figure  4l  . Load  Diagi’am  of  I Lng  Section. 
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The  average  pressure  acting  over  the  entire  wing  from 
assumption  one,  ist 

P.V.  = Ji. 

s 

where  S,  the  reference  area,  is  0.2393  ft^.  in  this  analysis, 
loads  were  considered  at  section  A-A,  where  the  bending 
moment  is  maximum,  and  at  section  B-B,  where  the  steel 
insert  ended  and  only  epoxy  carried  the  loads  as  shown  in 
Fig  42  and  Fig  43. 

Section  hzk 

Bending  Stress. 

For  both  the  forward  and  aft  swept  wings,  since 
the  planform  area  was  identical,  the  force  acting  outboard 
of  section  A-A  is  given  byi 

N*,.,  = (6) 

where  Sj^,  the  area  outboard  of  A-  A,  is  0.1786  ft^,  A 
normal  force  equal  to  64.2  Ibf  is  obtained  using  Eq  (6). 

For  the  bending  moment  at  section  A-A,  the  following 
equation  was  usedi 


where  Y from  Fig  42  and  Fig  43  is  2.57  inches. for  both 
wings.  This  results  in  a bonding  moment  of  165  inch-lbf 
at  this  section. 
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All  Dimonaiona  in  inchoa 

Ficure  42.  Schematic  of  Forward  Swept  Wing 
uaed  for  Streaa  Analysis 
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To  ovaluatG  the  strescotj  at  section  A~k,  the  steel 
insert  was  modeled  as  a thin  beam  shown  in  Fig  44 » 


An  expression  for  the  moment  of  inertia  about  the 
x-axis  is  given  byt 

I,  = 9^'  (7) 

\z 

where  a and  b are  the  dimensions  of  the  beam.  At  any  point 
on  the  cross-sectioni  the  bending  stress  0“  is  given  byt 

= Ki  Z 
Tx 

where  Z = b/2i  for  the  maximum  stress  at  the  surface. 
Substituting  Eq  (?)  into  Eq  (8)  yields  the  expression! 

(T  = rA  (9) 

a W" 

which  gives  the  maximum  stress  on  the  beam.  The  insert 
dimensions  at  section  A-A  are  3*125  by  0.125  inches. 

This  yields  a maximum  bending  stress  of  21,120  Ibf/in^ 
for  both  forward  and  aft  swept  wings. 


Shearing  Stress . 


The  maximum  shearing  stress i using  the  same 
thin  beam  model  in  Fig  44,  is  given  byi 


r = r 


f 3cl  t*  8 h \ (10) 

V ) 


where  T,  the  twisting  moment  isi 


T = ( X - 3.^)6  ) 


(11) 


for  each  wing  (Ref  llil94)i  Prom  Fig  42,  X = 5*18  inches 
for  the  aft  swept  and  results  in  a twisting  moment  of 
136  in-lbf  using  Eq  (11).  For  the  forward  swept  wing 
from  Fig  43,  X = 0.1  inch  and  results  in  T = 190  in-lbf. 
Using  Eq  (10) , a maximum  shearing  stress  of  8921  Ibf/in^ 
is  obtained  for  the  aft  swept  wing  and  12,464  Ibf/in^  for 
the  forward  swept  wing. 

Safety  Factors . 

For  the  insert  material, SAE  4l4o  steel, the  follow- 
ing properties  exist i 

PTy  (yield  tensile  stress)  =,  13^.000  psi 

('Ultimate  shearing  stress)  = 90,000  psi 

F^i^  (ultimate  tensile  stress)  = 150,000  psi 

Fg  (yield  shearing  stress)  “ 79,200  psi 

O' 

. (Ref  1 1I62) 

The  Factor  of  Gafoty  (FS),for  ultimate  and  yield, ic:  computed 
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using  the  following  interaction  formula  (Ref  5i20)i 

FS  * 1 (12) 


The  ultimate  PS  for  the  aft  swept  wing  was  5*81  and  5 >06 
for  the  forward  swept  wing.  The  yield  PS  was  5 >^3  for 
the  aft  wing  and  4,67  "the  forward  swept  wing. 

Section  B~B 

Bending  Stress . 

For  both  the  aft  and  forward  swept  wings*  the 
force  acting  outboard  of  section  B-B  is  given  asi 

^b-b  * • 

2 

where  the  area  outboard  of  B-B,  Sg  = 0.042?  ft  . This 
yields  a normal  force  of  15.4  Ibf  for  each  wing. 

For  the  bending  moment  at  section  B-B,  the  equation 

Ms-b  = l^e-s  y 

was  utilized  where  Y = 0.92  inches  for  both  wings.  This 
yields  a bending  moment  of  14.1?  in-lbf  for  each  wing. 

The  wing  section  at  B-B  is  also  modeled  as  a thin  beam 
with  dimensions,  a “ 3 *2  in.  and  b = 0.1655  in.  The  max- 
imum bending  stress  using  Eq  (9)  is  970  Ibf/in^. 

Shearing  Stress . 

For  both  wings,  the  twisting  moment  at  section 
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B-B  is  given  asi 

T = t X - i-es) 

Where  X = 2.6  in.  for  the  aft  wing  and  X = 0.80  in.  for 
the  forward  swept  wing.  The  twisting  moment  for  the  aft 
and  forward  swept  wings  is  11.55  in-lhf  and  15*6  in-lbfi 
respectively.  The  maximum  shearing  stress,  using  the 
Eq  (10)  and  the  dimensions  of  section  B-B,  was  407.6  Ibf/in^ 
for  the  aft  wing  and  550.3  Ibf/in^  for  the  forward  swept 
wing. 

Safety  Factors . 

The  aluminum  base  epoxy  used  to  mold  the  wings 
has  a bond  strength  of  4650  Ibf/in^.  Using  this  property 
and  Eq  (12),  the  ultimate  FS  at  section  B-B  for  the  aft 
swept  wing  was  computed  as  4.42.  For  the  forward  swept  wing 
the  FS  was  4.24. 

An  analysis  of  the  loads  on  the  fuselage  was  not 
considered  necessary  since  it  was  designed  for  loads  much 
greater  than  those  expected  in  this  test.  The  analysis  of 
this  component  can  be  found  in  Ref  5* 
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APPENDIX  B 

AERODYNAMIC  CHAftACTERISTICS  PRESENTED 
IN  GRAPHICAL  FORM 


The  results  of  the  wind  tunnel  tests  for  Mach  Numbers 
from  0.70  to  0.90  are  presented  graphically  in  this  appendix. 
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Ficuro  45 f vs  «<  for  M " 0.7 
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Figure  5^'  Cq  vs  c<  for  M = 0.9 
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Figure  55*  Cq  vs  Cj^  for  M = 0.7 
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Figure  62.  L/d  vs  for  M = 0.85 
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Figure  64.  l/D  vs  for  M = 0.9 
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Figure  65 « Cf/[  vs  for  M = 0.7 
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A FORWARD  SWEPT 


Figure  68.  vs  c<  for  M - 0.87 
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